By Mudit Sinhal 

For details see: Sinhal, M.: Quantum control of single molecular ions. PhD thesis, University of Basel (2021)

The file "StarkShift_allmixed_N2_HF.m" is the main file. Given any positive parity state in the X electronic state and a range of laser wavelengths it computes and can plot the ac-Stark shift and the scattering rate.
For example to calculate the ac-shift of the N = 0 and I = 2 nuclear spin isomer in the X^2\Sigma_g^+ state of N2+, run: "StarkShift_allmixed_N2_HF(0, 2);"  


The file "ac_Stark_allmixed_hfs.m" calculates the ac-Stark shift for a give transition, nuclear spin, hyperfine and projection of hyperfine quantum numbers. Here you can specify the laser power, laser polarization and laser focus as well. 


The file "honl_london_allmixed_hfs.m" calculates the honl london factor for positive parity X^2\Sigma_g^+ states to negative parity A^2\Pi_\Omega states in N2+. It is used for calculation of ac-shifts and scattering rates. 


The files "hamiltonian_XState_hfs.m" and "hamiltonian_AState_hfs_v2" calculate the transition energies and the mixing coefficients by solving the hamiltonian. The constants and hamiltonians are written in the file along with their references. 


The Files "Wigner3j.m" and "Wigner6j.m" are matlab functions for calculating the 3j and 6j functions. They are downloaded from matlab central
